2023年郑老师多门科研统计课程：多次直播，含孟德尔随机化方法

随着人口老龄化的进程，越来越多的研究关注如何在老年人中维持身心健康。其中，抑郁症是一种常见的情感障碍，与心理健康以及生理健康都有紧密的联系。虽然线性回归等传统的统计模型可以用于探究活动与抑郁症之间的关系，但并不能完全反映数据中可能存在的复杂结构。针对这个问题，聚类分析作为无监督式学习的重要工具，能够从数据层面上挖掘出样本之间的相似度和差异性，并将其分类为不同的群组。在本文中，我们将应用聚类分析方法，利用NHANES数据库中的数据，探究活动与抑郁症状以及衰老认知的关系，希望为后续老年精神保健实践提供一些启示性的结论。

2022年10月1日，一篇题为：Association of 24-Hour Activity Pattern PhenotypesWith Depression Symptoms and Cognitive Performance in Aging 的研究论文发表于《JAMA Psychiatry》（JAMA子刊）（IF=25.91）。

这项研究利用美国营养健康（NHANES）的数据，研究了在美国24小时活动模式表型与衰老抑郁症状和认知表现的关联。结果表明，美国大约五分之一的老年人可能被归类为活动模式较弱且安定较晚的亚组，大约十分之一的老年人可能被归类为活动模式较弱且活动持续时间较短的亚组。

一、摘要

背景：为了指导有针对性的机制研究和行为干预的开发，需要有关老年人24小时活动模式表型性质和普遍程度的证据，特别是其中与抑郁症状和认知有关的表型。本研究旨在确定具有相似24小时活动节律特征的老年人亚组，并描述相关的抑郁症状和认知表现。

方法：本次研究为横断面分析，使用了2011-2014年全国健康和营养调查(NHANES)中的加速度计研究数据。统计分析从2022年1月至3月。NHANES采用了多阶段随机抽样方式，旨在代表美国非机构化成年人口。主要分析对象为65岁及以上、具有加速度计和抑郁症状测量值的参与者，加权后代表约3200万老年人。主要的暴露变量是分组依据是运用潜在类别分析确定了子群体，这些子群体具有相似的24小时活动模式特征，这些特征是使用扩展余弦和非参数方法测量得出的。通过使用协变量调整样本加权回归，评估了子群体成员与以下两方面的相关性：（1）抑郁症状，定义为9项患者健康问卷(PHQ-9)得分大于等于10分；（2）至少存在心理测量上轻度认知障碍(p-MCI)，即综合认知表现得分低于平均值1个标准差。

结果：实际的聚类样本大小为1800(加权：平均[标准差]年龄72.9 [7.3]岁；57%女性参与者)。聚类确定了4个子群体：(1)677个早期起床/健壮型（37.6%）；(2)587个短期活动期/不易建模（32.6%）；(3)177个短期活动期/非常虚弱（9.8%）；(4)359个较晚结束活动/非常虚弱（20.0%）。PHQ-9得分达到或超过10分的患者比例在各组之间存在显著差异（聚类1：3.5%；聚类2：4.7%；聚类3：7.5%；聚类4：9.0%; χ2 P = .004）。至少存在p-MCI的患者比例在各组之间存在显著差异（聚类1：7.2%；聚类2：12.0%；聚类3：21.0%；聚类4：18.0%；χ2 P < .001）。其中有五个抑郁症状在子群组间存在显著差异。

结论：这项横断面研究的发现表明，在美国约5分之1的老年人可能被归为活动模式较弱且结束活动较晚的子群体，约10%的老年人可能被归为活动时间短且模式较弱的子群体。未来的研究需要探索与这些行为表型相关的生物学过程，包括为什么早期和强健的活动模式似乎是有保护作用的，以及改善异常的模式是否会改善结果。

二、研究结果

1. 活动模式的聚类识别（基线资料）

实际的聚类样本大小为1800(加权：平均[标准差]年龄72.9 [7.3]岁；57%女性参与者，43%男性参与者)。包括以下种族和族裔类别的参与者(按加权百分比计算)：亚洲人(聚类1：4.1%；聚类2：2.1%；聚类3：3.4%；聚类4：3.4%)，黑人(聚类1：5.3%；聚类2：6.6%；聚类3：12.4%；聚类4：13.5%)，西班牙裔 (聚类1：5.9%；聚类2：6.2%；聚类3：9.6%；聚类4：5.9%)，以及白人、其他人种或混合种族和民族 (聚类1：84.8%；聚类2：85.1%; 聚类3：74.6%；聚类4:77.3%)。参与者佩戴加速度计的平均(标准差)记录长度为6.9(0.2) 天 (范围4-7天)。BIC指出,最优模型具有4个子群体和相等形状之椭球协方差结构(见补充资料eTable 1和2)。我们将最大的子群体(677 [37.6%])标记为聚类1：早期起床/健壮型，因为他们平均上有：(1) 活动启动时间在早上7点之前（约早于平均值的三分之一SD） (2) 共15小时的活动时长；(3) 在所有指标中都体现出最强的活动模式。

第二大子群体(587 [32.6%])被标记为聚类2：短活动时期/不易建模，因为平均而言，他们与聚类1早期起床/健壮型群体有所区别，即开始日常活动时间较晚且结束时间较早(平均活跃周期持续时间=13.4个小时)，除虚F统计量外，他们的稳健性指标高于平均水平（其表明对24小时扩展余弦模型更劣）。最后的两个群体都表现出不规律的活动模式，多种指标均有所反映（伪F统计量、日际稳定性、日内变异性和相对幅度）。其中较小的那组占样本的9.8%（即1800人中的177人），被标记为聚类3：短活动时期/虚弱型，因为他们具有最短的平均活跃周期长度（11.8小时）。最后一组(359 [20.0%])被标记为聚类4：较晚结束活动/虚弱型，因为他们平均上具有最晚的估计活动偏移时间。表格1展示了各子群体的样本社会人口因素。

2.活动模式亚组中的抑郁症状分析

PHQ-9得分为10分或以上的患病率在不同组之间存在差异（聚类1，3.5％; 聚类2，4.7％; 聚类3，7.5％; 聚类4，9.0％；χ2 P = .004）（图2A）。总体而言，有5.4％的参与者PHQ-9得分为10分或更高。根据年龄、性别和种族调整后的数据，与聚类1早期起床/强壮型组相比，在聚类2活动时间较短/建模能力较差的组中出现PHQ-9得分为10分及以上的比例大致相同（比值比[OR]，1.43；95％ CI，0.76-2.67），而在聚类3活动时间较短/非常虚弱的组（OR，2.34；95％ CI，1.12-4.91）以及聚类4较晚结束活动/非常虚弱组（OR，2.91；95％ CI，1.51-5.60）中，此方面的比例均显著提高。

 此外，在9种个体抑郁症状中，各组之间也存在统计学上显著的差异（抑郁、睡眠障碍、疲劳/精力不足、内疚、自杀倾向）。与聚类1早期起床/强壮型组相比，有序回归和CI（图3）所指出的OR值表明（1）两个活动时间较短的组报告了更高的抑郁情绪水平（尽管这些组中最小的那个组的95％ CI包括1），（2）聚类4较晚结束活动/非常虚弱的组报告了更高的睡眠障碍水平，（3）所有三个群体都报告了更高的疲劳水平，（4）仅聚类3活动时间较短/非常虚弱的组报告较高的内疚水平，（5）两个活动模式较弱的群体（聚类3和4）报告了更高的自杀倾向水平。

3.活动模式亚组中的认知测试表现

需要注意，实际的未加权原始认知测试分数范围很广，表明样本涵盖了各种认知状态的参与者（见补充资料中的eTable 3）。至少有p-MCI的患病率在不同组之间存在显著差异（聚类1，7.2％; 聚类2，12.0％；聚类3，21.0％；聚类4，18.0％；χ2 P < .001）（图2B见上）。与聚类1早期起床/强壮型组相比，在所有3个组中（年龄，性别，种族，认知测验语言和教育已经调整），至少有p-MCI的患病率均较高（聚类2活动时间较短/建模能力较差的组：OR，1.84；95％ CI，1.06-3.20；聚类3活动时间较短/非常虚弱的组：OR，3.36；95％ CI，1.81-3.26；聚类4较晚结束活动/非常虚弱的组：OR，2.76；95％ CI，1.63-4.67）。

就连续认知表现结果而言，两组活动模式较弱的群体在所有三个测试上表现都不如其他群组（β系数> -0.25并且95％ CI不与0重叠）（图4）。相对于参照组，聚类2活动时间较短/建模能力较差的组仅在延迟记忆测试上表现较差。

https://mp.weixin.qq.com/s?__biz=MzAwOTYyMDY3OQ==&mid=2650405655&idx=4&sn=95f4717f40619275bf8172439fd2ac07&chksm=835186bfb4260fa9a25ee7eb3597d1e94aeae2322548bd000dae8be7fed4654559b0a0919de3&token=525144421&lang=zh_CN#rd

2023年统计服务

2023年统计服务开启！欢迎提前洽谈数据分析、科研合作服务